Dynamical ionization ignition of clusters in intense and short laser pulses 
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The electron dynamics of rare gas clusters in laser fields is investigated quantum mechanically 
by means of time-dependent density functional theory. The mechanism of early inner and outer 
ionization is revealed. The formation of an electron wave packet inside the cluster shortly after 
the first removal of a small amount of electron density is observed. By collisions with the cluster 
boundary the wave packet oscillation is driven into resonance with the laser field, hence leading 
to higher absorption of laser energy. Inner ionization is increased because the electric field of the 
bouncing electron wave packet adds up constructively to the laser field. The fastest electrons in the 
wave packet escape from the cluster as a whole so that outer ionization is increased as well. 



PACS numbers; 36.40.-c, 42.50.Hz, 33.80.-b, 31.15.-p 



I. INTRODUCTION 

Clusters bridge the gap between bulk material and 
single atoms. When placed into a laser field all atoms 
inside (not too big) clusters experience the same laser 
field, contrary to the bulk where a rapidly created sur- 
face plasma and the skin effect prevents the laser from 
penetrating deeper into the target. In rarefied gases, on 
the other hand, the laser can propagate but the density 
is too low to yield high absorption of laser energy and 
high abundances of fast particles. Hence, the absorption 
of laser energy is expected to be optimal in clusters. In 
fact, highly energetic electrons , ions Q, |^ Q , photons 

ly, , and neutrons originating from nuclear fusion f3] 
were observed in laser cluster interaction experiments. 
A prerequisite for clusters as an efficient source of en- 
ergetic particles and photons is the generation of high 
charge states inside the cluster. Several mechanisms for 
the increased ionization (as compared to single atoms in 
the same laser field) have been proposed. Boyer et al. 
suggested that collective electron motion (CEM) in- 
side the cluster is responsible for inner shell vacancies the 
corresponding radiation of which was observed in exper- 
iments 13 • Rose-Petruck et al. introduced the so- 
called "ionization ignition" (lONIG) where the combined 
field of the laser and the ions inside the cluster leads to 
increased inner ionization, i.e., electrons are more easily 
removed from their parent ion as if there was the laser 
field alone. The removal of electrons from the cluster as a 
whole is called outer ionization. An "outer-ionized" clus- 
ter will expand because of the Coulomb repulsion of the 
ions while a heated, quasi-neutral cluster expands ow- 
ing to the electron pressure. According to experiments 
the latter viewpoint of a nanoplasma (NP) ^ il2] seems 
to be appropriate for big clusters of about 10^ atoms or 
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more \13^ while numerical simulations indicate that the 
Coulomb explosion dynamics prevail for smaller clusters 
J5j- For recent reviews on the interaction of strong 
laser light with rare gas clusters see 0, ^3 ■ 

In this paper we are aiming at clarifying the early ion- 
ization dynamics in laser cluster interaction. Because 
the plasma has to be created first the NP model is not 
applicable. Although we shall find CEM to be indeed 
an important ingredient for the early ionization dynam- 
ics of clusters, the origin of the CEM should be explained 
rather than assumed beforehand. The IONIC model pre- 
dicts that "as the ionization proceeds beyond the first 
ionization stage, strong electric fields build up within the 
cluster that further enhance ionization" . Here we are 
interested in how these strong inner fields are generated. 
Once the cluster is charged after the first electrons are 
removed, the remaining bound electrons experience the 
attraction of neighboring ions. This attraction will be 
largest for bound electrons of ions sitting near the clus- 
ter boundary. The other ions will pull these electrons 
into the cluster interior. This force might be supported 
by the electric field of the laser, thus leading to further 
ionization. The IONIC model was successfully called on 
for interpreting the experimental results in [iSj. 

Although this IONIC mechanism is appealing, the de- 
tails of the ionization dynamics remain unclear. If the 
force exerted by the other ions is strong enough to ionize 
further it should be also strong enough to keep the al- 
ready freed electrons inside the cluster. Hence, it remains 
to be explained how the electrons are removed from the 
cluster as a whole (outer ionization). The electrons, af- 
ter inner ionization, may as well shield the space charge 
of the ions so that the IONIC mechanism would come 
to an end, and the NP model would take over even in 
small clusters. So, why does a strong electric field build 
up inside the cluster whose interplay with the electric 
field of the laser is constructive for inner as well as outer 
ionization? 

Classical particle methods are frequently applied to in- 
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vestigate the electron and io n dy namics of clusters in 
laser fields [H El El IM El El El since a full three- 
dimensional quantum treatment is out of reach with 
nowadays computers. In these classical accounts inner 
ionization was either accounted for by sampling the quan- 
tum mechanical electron distribution of the bound state 
by a microcanonical, classical ensemble of electrons, or 
by using ionization rates so that the electron dynamics 
was simulated only after inner ionization. Semi-classical 
Thomas-Fermi theory was employed in |15| to study the 
explosion dynamics of clusters consisting of up to 55 
atoms. In our work we use time-dependent density func- 
tional theory [3EHE3 since we are mainly interested in 
the early ionization dynamics of clusters were quantum 
mechanics is important. 

The article is organized as follows. In Section the 
numerical model is introduced. In Section [llll results are 
presented concerning the groundstate properties of the 
model cluster (|III All , the electron dynamics (jlll B|l . the 
formation of collective electron motion inside the cluster 
and outer ionization Ijlll C|) . the effect of the collective 
electron motion on inner ionization Ijlll D|l . the frequency 
dependence of outer ionization Ijlll and the influence 
of the ionic motion (|III F|l . Finally, we conclude in Sec- 
tion |IV| 

Atomic units (a.u.) are used throughout the paper 
unless noted otherwise. 



II. NUMERICAL MODEL 

Time-dependent density functional theory (TDDFT) is 
employed to study the ionization dynamics of small and 
medium size rare gas clusters in intense and short laser 
pulses. To that end the spin degenerate time-dependent 
Kohn-Sham equation (TDKSE) is solved. However, solv- 
ing the TDKSE in three spatial dimensions (3D) for 
rare gas clusters in laser fields is too demanding for to- 
days computers. One reason for this is that, contrary 
to metal clusters or fuUerenes [HI, the electrons in 
the ground state of rare gas clusters are not delocal- 
ized so that jellium models where the ionic background 
is smeared out and assumed spherical are not applica- 
ble. In order to make a numerical TDKS treatment fea- 
sible two simplifications were made. Firstjas in pre- 
vious studies of clusters in laser fields IH |^ Is^ . 
the dimensionality of the problem was restricted to ID, 
namely to the direction of the linearly polarized laser 
field described by the vector potential A{t) = A{t)e3. 
in dipole approximation ,33]. To that end the Coulomb 
interactions were replaced by soft-core Coulomb inter- 
actions, i.e., \r — r'\~^ [{x ~ x')'^ + «ee]~^^^ ^-nd 
|r - Rk\-^ [{x - Xkf + al^V'^/'^ for the electron- 
electron interaction and the electron-ion interaction, re- 
spectively. The smoothing parameters Oee and aei may 
be chosen to yield ionization energies similar to real 3D 
systems. The second simplification was the use of the 
exchange-only local density approximation (XLDA) so 



that V^^[n{x,t)] = Vyii^YyA[n{x,t)\ = -a(3n(a;,t)/7r)i/3 
where n{x, t) is the electron density. The pre-factor a 
would be unity in full 3D XLDA calculations. In the ID 
model it may be chosen to yield satisfactory ground state 
properties (we used a = 3/4). The ionic potential was 

Mon(a;) - - J2k=i Z[{x - Xkf + al,]-^/'^ with constant 
nearest-neighbor distances Xk+i — Xk — d ^ 2rws where 
r-ws is the Wigner-Seitz radius. One may look at the ID 
ion chain as representing those A^ion ions of a 3D spheri- 
cal cluster which are situated along the diameter parallel 
to the linearly polarized laser field. The cluster radius 
then is i? w (A^ion — 1)^/2, and the number of ions in the 
real 3D cluster Nsd = R^/^ws would be « (Mon - 1)^- 
The TDKSE for the A^ks Kohn-Sham (KS) orbitals 
'i'jix, t), j = 1,..., Nks reads 

i^*j(a^,0 =^KS*j(x,0 (1) 

where 

i?KS = -^V2+C/+VScLDA + Mon+Mascr. (2) 

From the doubly spin-degenerate KS orbitals the total 
electron density n{x, 

t) = ^Ej=^\^]ix,t)\^ is calcu- 
lated. U{x,t)^ /da;'n(a;',t)[(a;-a;')2-f-a2j-i/2 isthclD 
Hartree potential (accounting for the electron-electron 
repulsion), and Mascr = —^A{t)dx governs the interaction 
with the laser (taken in velocity gauge with the purely 
time-dependent term ~ transformed away). Eq. 
was solved using the Crank-Nicolson method with the ac- 
curacy of the spatial derivatives boosted to fourth order. 
Since the Hamiltonian Q itself depends (through the 
density) on the KS orbitals a predictor-corrector method 
should be used in combination with the usual Crank- 
Nicolson time-propagation. In practice, however, the 
evaluation of the Hamiltonian |(2Jl using the density from 
the previous time step is usually accurate enough. 

For the ion-ion interaction a soft-core potential 
Z'^[{Xj — Xk)'^ + a|]^^/^ was assumed as well. The ions 
j — 1, . . • , A'ion were moved according to their classical, 
non-relativistic equations of motion, 

N- 

ion Y — Y 

_ f n{x,t)iXj-x) 

using the Verlet algorithm. 

The results which will be presented in the Sections 
nil Al ITTTeI were obtained for an ion mass M — 131 • 1836 
(Xe atom) and an = 1. Due to the short laser pulse 
durations (< 40 fs) and the modest charge states < 4 
created, the ionic motion did not appreciably affect the 
ionization dynamics of our ID model. However, since in a 
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real 3D cluster consisting of (iVion — 1)^ ions the Coulomb 
explosion would be more violent than in our ID model 
with A^ion ions (of the same ion mass and charge state) 
the issue of ionic motion will be discussed separately in 
Section ImFl 

Finally, before presenting our results, we briefly relate 
our approach to the work in Refs. I29I |3 lL l3 ^ . Instead 
of the simple XLDA, Veniard et al. used in their 
ID TDDFT approach the more advanced exchange po- 
tential proposed by Krieger, Li, and lafrate (KLI) |3J| 
but restricted their studies to ion chains up to 7 fixed 
atoms with one active shell only. Besides the extensive 
study of the harmonic radiation emitted by the ID clus- 
ter, Veniard et al. focussed on the dependence of ion- 
ization on the internuclear distance and on the number 
of ions in the chain. This is different from our goal to 
highlight the dynamics of lONIG. The methodological 
approach instead is very similar. The only differences are 
the mobile ions and the choice of simple XLDA instead 
of the KLI exchange potential. The latter simplification 
gave us the freedom to simulate bigger clusters and more 
active shells per atom. 

The ID models in and 0| were called on to study 
the explosion dynamics of Xe clusters. Since for this 
purpose a full TDDFT treatment with many KS orbitals 
is too demanding even in ID, further simplifications had 
to be adopted. In |33| a single orbital, representing all 
electrons, was introduced while in (2^ time-dependent 
Thomas-Fermi theory was used. Our main emphasis in 
the present paper is complementary to this work since we 
are not dealing so much with the Coulomb explosion of 
the cluster but are predominantly interested in the early 
electron dynamics. 



III. RESULTS AND DISCUSSION 

A. Groundstate properties 

Let us consider a chain of Mon = 9 ions with nearest- 
neighbor distance d = 8 and nuclear charge per ion 
Z = A. Hence, the ground state of the neutral cluster 
consists of = 36 electrons. The smoothing parameters 
for the soft-core Coulomb interaction were simply chosen 
to be Oee = flei = flji ~ 1- In Fig. n the ground state 
density and the various contributions to the total po- 
tential are plotted: the ion potential VJoid the Hartree 
potential U, and the exchange potential Vxlda- Al- 
though the ionic potential alone has its absolute mini- 
mum at the central ion the total potential, including the 
classical Hartree-repulsion and XLDA, consists of nine 
almost identical and equidistantly separated potential 
wells, each locally similar to that of the corresponding 
individual atom. Consequently, the ground state density 
displays nine almost equal, well localized peaks, and the 
energy levels of the cluster are approximately at the same 
positions as those of the atom, namely around —1.17 for 
the 2A'ion — 18 inner electrons, and around —0.23 for 



FIG. 1: Ground state electron density (orange, dashed- 
dotted) and the various contributions to the total potential 
(red, solid): ionic potential Vion (black, solid), Hartree po- 
tential U (blue, dotted), and XLDA potential VxLDA (green, 
dashed). 



FIG. 2: Logarithmic density log n{x, t) vs. space and time 
for (a) the cluster with Ni^^^ = 9, Z — 4 (b) the single atom, 
and (c) a cluster made of noninteracting atoms, all with an 
electron dynamics as shown in (b). The laser parameters were 
uj = 0.057, E = 0.033, (3, 8, 3)-pulse. 



the 18 outer electrons (note, that in the ID model there 
are only 2 electrons per shell). Since in XLDA Koop- 
man's theorem is usually not well fulfilled the energy of 
the highest occupied orbital does not equal the ioniza- 
tion energy for removing the outermost electron. Cal- 
culating the ionization energy from the difference of the 
total energy of the neutral and the singly ionized cluster 
Afciuster = 0.26 was obtained. Doing the same for the 
single atom led to A^atom = 0.49 which is a reasonable 
value for rare gas atoms. 



B. Electron dynamics 

Let us start by comparing the electron motion in the 
cluster with that of a single atom. In Fig. [21 contour plots 
of the logarithmic density are shown vs. space and time. 
The results for the full cluster (a), the single atom (b), 
and an artificial cluster made of noninteracting atoms (c) 
in a laser pulse of frequency oj — 0.057 (A = 800 nm) and 
rather modest field amplitude E ~ 0.033, corresponding 
to « 3.9 X lO^^Wcm^^, are shown. The laser field was 
ramped up linearly over 3 cycles, held 8 cycles constant, 
and ramped down again over 3 cycles (hereafter called 
a (3, 8, 3)-pulse) so that the pulse duration was w 37 fs. 
The density in the contour plot (c) was calculated by 
assuming that all of the A^ion = 9 atoms behave as the 
single atom in plot (b). 

It is seen that the electron dynamics of the cluster 
(a) and the noninteracting atoms (c) differ significantly 
from each other already at t « 150 because more electron 
density leaves to the right at that time instant in (a) than 
it does in (c) . This behavior of stronger ionization of the 
cluster than of the individual atoms continues during the 
subsequent half laser cycles. 

A qualitatively different electron dynamics inside the 
cluster emerges for t > 400. While in the cluster (a) 
an accumulation of electron density bounces from one 
boundary of the cluster to the other such an electron 
dynamics, of course, cannot build up in the ensemble 
of independent atoms (c). The formation of an electron 
wave packet which travels through the entire cluster in 
step with the laser field is remarkable in view of the fact 
that the excursion of a free electron in the laser field 
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E = 0.033 amounts to x = 10.2 only while the diameter 
of the cluster is 2R « 64. By varying the cluster size 
(results for A'^ion = 17 will be presented in the following 
subsection) it was found that the formation of the bounc- 
ing wave packet at laser intensities where E/cu'^ < R/2 
is a robust phenomenon, not sensitive to the cluster pa- 
rameters. Test runs with more advanced exchange po- 
tentials such as the Slater potential with self-interaction 
corrected XLDA and the KLI potential were per- 
formed to ensure that the coherent electron motion is 
not an artifact of plain XLDA. However, the wave packet 
formation is sensitive to the laser frequency, as will be 
shown in subsection IIII El 



C. Formation of collective electron motion inside 
the cluster and outer ionization 

It is useful to study the phase relation between the os- 
cillating electron density inside the cluster and the laser 
field in order to understand the formation of the elec- 
tron wave packet with unexpected large excursion ampli- 
tude. Given a laser field E(t) of optical (or lower) fre- 
quency, the polarization of an atom is ^ —E{t) because 
the bound electrons are able to follow adiabatically the 
force exerted by the field. Hence, the phase lag of the po- 
larization with respect to the laser field is A(/) = tt. Free 
electrons, on the other hand, oscillate perfectly in phase 
~ E{t) so that A(/) =: 0. Energy absorption from the laser 
field is low in both cases because / dtx{t) ■ E{t) k, 
(with x(t) the expectation value for the position of an 
electron) . During ionization there is necessarily a transi- 
tion from A(/() = TT to where energy absorption can take 
place. During the ionization of atoms this transition oc- 
curs rapidly while in clusters, after inner ionization, the 
electrons may be still bound with respect to the cluster 
as a whole. Hence, the free motion of electrons ~ E{t) 
comes to an end at latest when they arrive at the cluster 
boundary. There, they either escape from the cluster, 
contributing to outer ionization, or they are reflected so 
that their phase relation with the driving laser is affected, 
leading on average to an enhanced absorption of laser 
energy. Although collisions of the electrons with ions 
are included in our TDDFT treatment the effect of the 
boundary on the electron dynamics clearly dominates. 
The unimportance of electron ion collisions in medium 
size and small clusters was pointed out in |14 | while the 
relevance of boundary effects was recently afltirmed in 
within the NP model "s^. 

In Fig. 12^- it is seen that some electrons enter about 
ten atomic units into the vacuum before they are pulled 
back by the cluster charge. This is reminiscent of what 
in laser plasma physics is called Brunei effect |3S|, "vac- 
uum hea ting " p39| , or, more expressively, "interface phase 
mixing" jlj. Thanks to the fact that the fast electrons 
leave the cluster (and slow electrons are accelerated) a fil- 
ter effect comes into play so that a wave packet can form 
that oscillates with the laser frequency and an excursion 



FIG. 3: Dipole 

^^inner and number of electrons inside the clus- 
ter A'^innoi. vs. time for different laser and cluster parameters. 

(a) E = 0.033, uj = 0.057, iVion = 9, Z = 4 (as in Fig. 

(b) E = 0.114, UJ = 0.057, iVion = 9, = 4; (c) i = 0.114, 
u = 0.057, iVion = 17, Z = 2; (d) E = 0.099, u = 0.18, 
A^ion = 9, Z = 4. The course of -E{t) (a (3, 8, 3)-pulse) is in- 
cluded in the a;innor-plots (dotted in green) for distinguishing 
motion in phase with —E{t) (green bar at the top edge of the 
panel) and motion approximately 7r/2 and tt out of phase (red 
and blue bar, respectively). For comparison, A^jon times the 
result for the single atoms are included in the Ajnner-plots. 



amplitude of about the cluster radius R. 

In order to underpin this scenario the dipole of 
the electron density inside the cluster XinnGr(i) ~ 
J^^^'^^^.^dx xn{x,t) was calculated for several laser and 
cluster parameters. The results are shown in Fig. |3| In 
panel (a) the laser and cluster parameters were the same 
as in Fig. [5] One sees that during the first few laser pe- 
riods the electrons indeed move ~ —E{t) (green, dotted 
curve) as indicated by the first, green bar at the top of 
the Xinner-plot. Then, the electrons inside the cluster get 
out of phase with the laser for about nine cycles (red bar) 
so that A0 K, Ti/2. During this period the dipole ampli- 
tude Xinnor IS particularly high |4j| , and the wave packet 
bouncing inside the cluster is clearly visible in Fig. |21 Fi- 
nally, towards the end of the laser pulse the phase relation 
of the few electrons which were removed from their par- 
ent ions but did not make it to leave the cluster becomes 
that of free electrons (blue bar), i.e., Xinnci{t) ^ E{t). In 
the lower plot the number of electrons inside the cluster 
A^innor IS plotted and compared with A^ion times the result 
for the single atom. It is seen that during the first phase 
(green bar) ionization of the cluster proceeds similar to 
the single atom case. However, when the phase lag is 
shifted to A0 = 7r/2 the cluster continues to ionize while 
the single atom ionization comes to an end. 

In Fig. I^ls the same is shown for a higher laser inten- 
sity. Essentially, the phase A0 behaves in the same way 
but this time, due to the stronger laser field, ionization 
of the outer shell is almost completed during the first few 
laser cycles. Hence, the final average charge state is al- 
most the same for the cluster and the single atom. The 
period where the phase lag is about A0 « 7r/2 lasts only 
a few laser cycles (red bar) and so does the wave packet 
motion inside the cluster. In Fig. O; the result for a big- 
ger cluster (A^ion = 17, Z = 2) is presented, revealing a 
qualitatively similar scenario as in (a) with the bounc- 
ing wave packet surviving for about 9 cycles. Finally, in 
panel (d) a higher laser frequency was used (w = 0.18) 
while keeping the cluster parameters as in (a) and (b). 
The A^innor-plot rcvcals that the single atom ionizes more 
efficiently than the cluster for these laser parameters. We 
will come back to the frequency dependence of outer ion- 
ization in subsection IIII El 
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D. Effect of the collective electron motion on inner 
ionization: dynamical ionization ignition 

The formation of collective electron dynamics as ex- 
posed in the previous subsection explains how the ab- 
sorption of laser energy is increased due to a phase shift 
into resonance with the driving field, and how the elec- 
trons, after inner ionization, are efficiently transported 
out of the cluster (outer ionization). The increased in- 
ner ionization still remains to be analyzed. lONIG states 
that the presence of the other ions is responsible for the 
more efficient removal of bound electrons. This is because 
two neighboring ions form a potential barrier (cf. the Vion 
curve in Fig. ^ through which an electron may tunnel 
when the whole cluster is submitted to an electric field 
so that the entire potential is tilted. However, we found 
from our numerical studies that this energetic advantage 
of a bound electron inside the cluster as compared with 
an electron in the corresponding single ion is not very 
pronounced for medium size and small clusters at mod- 
erate charge states. Instead, we propose a dynamical 
version of lONIG where the previously introduced collec- 
tive electron motion plays an important role. Coherent 
electron motion was suggested to be responsible for inner 
shell vacancies in Xe clusters, leading to x-ray emission 

^3 ■ III oui' numerical model there are only two shells 
and we did not find particularly high line emission from 
the cluster as compared to the single atom. However, it 
is possible that, owing to the lack of dynamical corre- 
lation in XLDA, the interaction of the electrons in the 
wave packet with the still bound electrons is underesti- 
mated in our model. Thus, dynamical lONIG might be 
the mechanism behind the experimental results reported 
in 0. 

Despite the fact that in our mean field approach there 
are no "hard" collisions of the electrons in the wave 
packet with the still bound electrons, the electric field 
associated with the oscillating electron packet already 
enhances inner ionization. In Fig. 0] two snapshots of 
the total effective potential and the electron density are 
presented for the same laser parameters as in Fig. |5| In 
panel (a), electron density and total potential are shown 
for a time where the electron wave packet is close to the 
left cluster boundary (the red bar at the bottom of the 
density plot indicates x-mnorit))- The contour plot (b) 
shows for each ion at position Xi the "difference density" 

n(a;, t) — Ni^ner/Nion which indicates whether 
there is a lack of electron density at that position inside 
the cluster (black and dark colors) or whether there is 
excess density (yellow and light colors) compared to the 
average density Mnnor/Mon- 

In panel (a), at time t = 662.5 the wave packet is close 
to the left boundary of the cluster while there is a lack 
of electron density near the right boundary. This charge 
distribution leads to a force Fjnt on the other electrons 
(pointing to the right) and therefore increases the ion- 
ization probability. The electric field of the laser instead 
is close to zero so that i^iascr is small. The situation ap- 



FIG. 4: Snapshots of the electron density and the total po- 
tential (plus the laser electric field potential alone) at time 
(a) t = 662.5 and (c) t — 697.5. The bar at the bottom of 
the density plot indicates a;inncr(i)- Arrows, and — in the 
potential plots illustrate the forces Fiag^j., Fj^t exerted by the 
laser field and the space charge, respectively. The contour 
plots (b) and (c) show whether excess (light colors) or lack 
(dark colors) of electron density (with respect to the average 
density A'^inncrZ-^ion) prevails at the position of an ion. The 
black horizontal lines indicate the times where the snapshots 
were taken. 



FIG. 5: Z = 4 minus the electron density, integrated ±d/2 
around the central ion in the cluster to determine the number 
of removed electrons as a function of time (solid curve, red). 
The result for the single atom is also shown (dotted, black). 

proximately a quarter of a laser cycle later is shown in 
panel (c). The wave packet is at the center of the cluster, 
moving with maximum velocity to the right and repelling 
bound or slow electrons in front of it. The force i^iasor is 
close to its maximum value at that time, pointing into 
the direction in which the wave packet moves. The ion- 
ization probability is, again, greater than with the laser 
field alone. Thus, during the course of a laser cycle the 
total force i^iascr + ^int clcarly leads to higher ionization 
as if there was the laser field only. 

The fact that the electron wave packet dynamics in- 
deed increases inner ionization is underpinned by Fig. [S] 
where the number of electrons in the region ±d/2 around 
the central ion in the cluster was subtracted from the ini- 
tial value Z = A and is compared with the corresponding 
single atom result. The laser parameters were the same 
as in Figs.[21and^ While for the single atom the ioniza- 
tion is completed for t > 700 the average electron density 
around the central ion in the cluster is still decreasing. 
When the electron wave packet sweeps over the central 
ion the density is temporarily increased, leading to local 
minima in the curve of Fig. |S1 When the wave packet 
is closest to one of the two cluster boundaries the lack 
of electrons around the central ion is maximal. The ab- 
solute increase of this maximum each half cycle means 
ongoing inner ionization. In contrast, the single atom, 
where only the laser field is present but no wave packet 
can form, does not ionize any further. 

E. Dependence of outer ionization on the laser 
frequency 

The interaction of the model cluster with Nion = 9 was 
investigated for the two different laser frequencies coi = 
0.057 and ujh = 0.18 (corresponding to 800 and 254 nm, 
respectively) and laser intensities between 4 x 10^^ and 
lO^^ Wcm-2. The pulses were of (3, 8, 3)-shape for both 
frequencies, that is, the pulse durations were T; « 37 fs 
and Th = 12 fs, respectively. After the laser pulse, the 
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FIG. 6: Average charge state vs. laser intensity of the clus- 
ter (solid lines) and the individual atom (dotted) for the two 
different frequencies uii = 0.057 (If, drawn red) and LUh = 0.18 
(hf, drawn blue). See text for discussion. 



average charge state in the cluster Zav = Z — iVinncr/Mon 
was calculated. Note that Zav only yields information 
about outer ionization for it does not distinguish between 
electrons that are still bound to their parent ions and 
those which move inside the cluster. 

In Fig.Elthe average charge state is plotted vs. the laser 
intensity for the two frequencies Uh and lui. The results 
for the single atom are also shown. As discussed in the 
previous subsections, it is seen that in the low frequency 
case the atoms in the cluster are stronger ionized than 
an individual atom in the same laser field. Both charge 
states come close only for Zav = 2, that is when the 
two electrons of the first shell are removed but the two 
electrons of the next shell are still strongly bound. The 
stepwise increase of the charge state due to the electronic 
shell structure is very pronounced in the low-frequency 
cluster case as well as for the single atoms at both, low 
and high frequency. The cluster in the high-frequency 
field instead shows a very different behavior: between 
lO^'^Wcm"^ and the threshold to the inner shell at w 
5 X lO-'^^Wcm^^, the charge state of the single atom is 
higher than the average charge state in the cluster. 

In Fig.|3the dynamics of the cluster electrons is shown 
for a (3, 8, 3)-pulse of frequency tUh = 0.18 and peak field 
amplitude E — 0.099, corresponding to an intensity of 
3.44 X 10-^'*Wcm~^. It is seen from Fig. that for this 
intensity the single atom ionizes more efficiently than the 
cluster as a whole, contrary to what happens at the lower 
frequency coi = 0.057. Fig. |7| reveals that the electrons, 
although removed from their parent ions, mostly remain 
inside the cluster. From plot (c) one infers that if the 
atoms inside the cluster were independent there would be 
a strong electron emission for 100 < t < 350. The emit- 
ted electrons have sufficient high kinetic energy to escape 
from their parent ion (and, thus, from the "independent 
atom" -cluster). Contour plot (a), instead, shows that in 
the real cluster a significant fraction of the electrons near 
the cluster boundaries return due to the space charge 
created by all the ions. A wave packet dynamics that 
could enhance outer ionization, as in the low frequency 
result of Fig. El does not form. We attribute this to the 
fact that the initial inner ionization occurs less adiabatic 
(multiphoton instead of tunneling ionization). Moreover, 
the excursion x — 3.06 < d/2 is too small to trigger 
any collective motion. The dynamics inside the cluster is 
rather "splash-like," as can be inferred from the strongly 
fiuctuating electron density between the ions in Fig. [7^. 
Hence, although at high laser frequencies inner ionization 
is high, outer ionization remains low since there is not 
the coherent electron dynamics supporting outer ioniza- 
tion. Consequently, for creating quasi-neutral nanoplas- 



FIG. 7: Same as in Fig.|5|but for the higher frequency uih ~ 
0.18 and E = 0.099. Logarithmic density logn{x,t) vs. space 
and time for (a) the cluster, (b) the single atom, and (c) a 
cluster made of noninteracting atoms. 

mas and suppressing Coulomb explosion the use of high 
frequency lasers is favorable. Reduced ionization of clus- 
ters in laser fields of, however, many times higher fre- 
quency (to be generated by x-ray free electron lasers in 
the near future) was also found in the numerical simula- 
tions of [23|. On the other hand, in the soft-x-ray FEL 
experiment performed by Wabnitz et at |42i] increased 
ionization of Xe clusters (as compared to single atoms) 
was observed at 98 nm wavelength, 100 fs pulse duration, 
and intensities up to 7 x IQ-'^^Wcm"^. The reason for 
this unexpected behavior at short wavelengths is not yet 
clear. 

For the high frequency ujh and laser intensities / > 
3 X 10^^ Wcm~^ the average charge state in the cluster 
overtakes the charge state of the single atom (see Fig.EJ. 
It was therefore interesting to check whether under those 
conditions a wave packet also forms at the higher fre- 
quency LUh- This is indeed the case. However, due to the 
rapid ionization the wave packet dynamics lasts only a 
few laser cycles (or less). 



F. Influence of ionic motion on the wave packet 
formation 

So far the ion mass was set to M = 131 • 1836 (Xe 
atom). This high mass lead to no appreciable ionic mo- 
tion in the ID cluster during the pulse. Even in the 
"worst case" where all active electrons are removed at 
t = so that the ionic charge Z = A is not screened at 
all, the iVion = 9-cluster radius increases only from 32 a.u. 
to 33.8 a.u. within the 37 fs pulse duration. However, the 
corresponding spherical 3D cluster with (iVion — 1)'^ Xe- 
ions of charge state Z — -i doubles its initial radius within 
the same time |4^. Clearly, the Coulomb explosion is 
underestimated in the ID model because in 3D an ion 
sitting at the cluster surface "sees" not only the charges 
of the ions aligned along the laser field polarization direc- 
tion but essentially a charged sphere containing all the 
other (A^ion -1)^-1 ions. 

Hence, in order to allow the ID ion chain to explode 
similarly to the corresponding 3D cluster one has to re- 
duce the ion mass. To obtain the results shown in Fig.|S| 
M = 1836 was set. In the "worst case" introduced in 
the previous paragraph the radius of the ID cluster now 
expands from 32 to « 130 a.u. within 37 fs. However, 
in the early stage when the wave packet formation in- 
side the cluster takes place the charge states of the ions 
are still low and the cluster radius remains close to its 
original value. This is clearly visible in Fig. where 
between t = 300 and t = 800 the bouncing wave packet 
can be easily identified. In Fig. |Sl3,c the dipole Xinncr 
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FIG. 8: Result for E = 0.062, uj = 0.057, iVion = 9, Z = 4, 
(3,8,3)-pulse, and mobile ions (M = 1836). (a) Logarithmi- 
cally scaled electron density, (b) a;inncri said (c) A^inncr (a-s in 
Fig. 1^ . The blue and broken curves in (b) and (c) are the 
results for immobile ions. 



develop in such intense fields the stronger ionization of 
the cluster still relies on the fact that the electron density 
inside the cluster moves temporarily with the appropriate 
phase lag necessary for efficient absorption. 



FIG. 9: The same as in Fig. |Hlbut for E = 0.57. 

and the number of electrons inside the cluster A^inncr are 
presented, as in Fig. |51 For comparison, the results for a 
computer run with immobile ions is included (blue and 
broken curves). The differences are small. Only the ion- 
ization at the end of the laser pulse is slightly reduced in 
the case of mobile ions. This is expected because the ions 
gain their kinetic energy at the expense of the electrons. 

At higher laser intensities the higher charge states are 
generated earlier. In these cases the Coulomb explosion is 
more violent but the wave packet mechanism contributes 
little to the total ionization of the cluster anyway. In 
Fig. El the ionization and Coulomb explosion scenario for 
the highest laser intensity / — 1.14 x lO^^Wcm"^ con- 
sidered in this paper is presented for the same cluster 
parameters as in Fig. |H1 The outer shell of all atoms is 
depleted rapidly already during the first laser cycle both 
for the cluster and the isolated atom. On average three 
quarters of the electrons in the next shell are depopulated 
during the remainder of the laser pulse in the case of the 
cluster with mobile ions. The ionization probability of 
the isolated atom is slightly less. The ionization of the 
cluster with heavier ions is, again, higher. 

Due to the rapidly increasing charge states of the 
ions no formation of an electron wave packet moving 
Ac/) — 7t/2 out of phase with respect to the driving laser 
field for several laser cycles can be inferred in Fig.|^. In- 
stead, the typical motion ~ E{t) of free electron density 
in the laser field E{t) seems to dominate the dynamics. It 
is interesting to observe that the removal of the electrons 
in the first shell is even more efficient for the isolated 
atom than in the cluster (see Fig. ^ around t — 180). 
However, shortly after "cracking" the second shell the 
cluster ionization overtakes the single atom result. At 
that time the Xinnor-plot in Fig. |^ reveals two extrema 
(marked by arrows) which are approximately A(f) — tt/2 
out of phase with the driving field. Thus, although a dc- 
phased wave packet dynamics over several cycles cannot 



IV. CONCLUSION 

The ionization dynamics of a one-dimensional rare gas 
cluster model in intense and short laser pulses was inves- 
tigated by means of time-dependent density functional 
theory. An electron wave packet dynamics was found 
to build up inside the cluster when the laser intensity 
/ = was sufficiently high for modest inner ionization 
but not that high that all electrons of an atomic shell are 
freed within a few laser cycles. The electron wave packet 
is driven into resonance with the laser field through the 
collisions with the cluster boundary. The phase lag be- 
tween the bouncing electron wave packet and the laser 
field then is it/ 2 so that the absorption of laser energy 
is particularly high. The fastest electrons in the wave 
packet escape from the cluster (outer ionization). The 
electric field of the bouncing electron wave packet adds 
up constructively to the laser field, thus enhancing inner 
ionization. This effect was called dynamical ionization 
ignition. It is a robust phenomenon with respect to the 
cluster size and, since it occurs during the early ioniza- 
tion stage, it is not affected by ionic motion. However, 
with increasing laser frequency (keeping the laser inten- 
sity fixed) the mechanism is less efficient. 

We expect the wave packet scenario being valid also 
in real, three-dimensional rare gas clusters. Due to 
the spherical cluster-vacuum boundary the wave packet 
should assume a sickle-like shape in that case. However, 
in order to verify this, studies with higher-dimensional 
cluster models will be pursued in the future. 
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